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Abstract. The problem of determination of the velocity anisotropy pa- 
rameter from star density profiles in globular clusters is considered. The 
calculations are performed for 10 globular clusters. The simple way for 
study of the dependence of anisotropy parameter on radial distance is 
used. 



1. Introduction 

The anisotropic models occup an important place among various realistic mod- 
els for study of the dynamical evolution of globular clusters. In frame of this 
models it have been work out the many evolution theories of the kinematics 
and dynamics of spherical self-gravitating systems (e.g., Bettwieser &: Spurzem 
1986; Spurzem 1996; Louis 1990; Solanes & Salvador-Sole 1990; Takahashi 1996). 
At the same time have been accumlated many observational data. These data 
mostly concerned star counts and surface brightness profiles, proper motions, 
stellar radial velocity dispersion, etc. (Meylan & Heggie 1997). Despite of the 
presence such great quantity information, still there are difficulties in compari- 
sion problem the observational data with results of numerical simulations in the 
frame work of anisotropic models. One of causes of this phenomenon perhaps 
is concluded in that for their performe we must have at least the value of the 
anisotropy parameter A = 2 — 2of /of , where of and of are radial and tangential 
components of velocity dispersion respectively (0 < A < 2). In general case it 
depends on radial distance and time. The aim of this work is to determine the 
value of anisotropy parameter for 10 globular clusters, to try find out the its 
dependence on limiting magnitude and on radial distance. 



2. The Method 

We have used the method that was proposed by Agekian and Petrovskaya (1962). 
Let us assume: 1) the star ensemble has an equal-mass spectrum; 2) the cluster 
is spherically symmetric and is in stationary state in the regular force-field; 3) 
the cluster consists of single stars; 4) of = const and of = const along radius 
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(so A = const). Then the cluster is described by means of the equation 

where £ and Z) are undimensional radial distance and undimensional spatial 
density. A theoretical surface mass density profile is determined as 

oo 

FtheoAr) = -2(3 J ^?-(r/a) 2 D' (£) # - F , (2) 

r/a 

where a and (3 are scale factors, Fq is constant taking into account the infinity 
radius of the model and D' (£) is the solution of (1). From assumption 1) it 
follows that p oc n, where n is spatial number star density. Then Ftheor{ r ) is 
theoretical surface number density of stars. Here constants a, (3 and Fq are 

N 

determined by minimization of the function $>(A) = J2 Wtheor{ r i) — Fobs( r i)\ 

i=i 

where N is number of star counts zones and F fy S (r) is observed surface number 
density of stars. That value of A that gives minimum $(A) is adopted as most 
probable value of A. 



3. Results and Discussion 

We determined the value of anisotropy parameter A for following 10 globular 
clusters: M3, M12, M13, M15, M55, M56, Palomar 3, Palomar 14, NGC5139 (u 
Cen) and NGC6535. The results are given in Table 1. For most clusters we have 
A = 0, i.e. they have almost isotropic velocity distribution. It is interesting that 
for all clusters theoretical number density profiles are in good agreement with 
observational ones. But the theoretical density profile in the halo almost in all 
cases is lower than observational one. Perhaps the cause is in the assumption 
about equal-mass. The fact is, as a result of mass segregation the low-mass 
stars prevail in the halo and consequently we overestimated mass density in this 
region. In reality the mass density less than we suppose. 

For M56 at limiting magnitude B = 21™5 we found A = 0.8. Hence in 
this cluster the velocity distribution is strongly anisotropic yet. Besides in order 
to find out an influence of limiting magnitude to value of A, we calculated it 
at B = 19^0 and found A = 0.3, i.e. if in this case the limiting magnitude 
decreases, then the anisotropy parameter decreases also. Perhaps it is because, 
in first case {B = 21™5) it is expected the influence of low-mass stars to value 
of A that generally are in halo and move along elongated orbits. 

In order to find out the possible dependence of the anisotropy parameter on 
the radial distance we performed for NGC5139 (to Cen) additional calculations 
in three separate regions. The results show that in the central region A = 0, in 
middle region A = 0.1 and in the halo A = 0.3. In this case the agreement of 
the theoretical and observational density profiles essentially was improved than 
one for whole cluster. 

We think this way allowes to solve partly the problem of dependence of 
the anisotropy parameter on radial distance for any globular cluster. If we will 
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Table 1. Value of anisotropy parameter for 10 globular clusters. 



Cluster 


A 


Limit, magnit. 
or Exp. time 


Source of 
observational data 


Note 


M3 





10 min 


King et al. 1968 
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9 min 


King et al. lybo 
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10 min 


King et al. 1968 
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30 min 


iVing et al. lybo 
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2 min 


King et al. lybo 
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30 min 


Ivmg et al. lybo 
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Harris & 

van Hpti Rprcr 1 QR4- 




NGC6535 





B=20 m 


Peykov & 
Roussev 1988 




M56 


0.8 


B=21 m 5 


Peykov & 
Roussev 1986 




-II- 


0.3 


B=19 m 


Peykov & 
Roussev 1986 




NGC5139 





30 min 


King et al. 1968 


6.27 < r' < 23.82 


-II- 





-II- 


King et al. 1968 


6.27 < r' < 10.22 


-II- 


0.1 


-II- 


King et al. 1968 


10.22 < r' < 14.75 


-II- 


0.3 


-II- 


King et al. 1968 


14.75 < r' < 23.82 



make use of the more precise observational data, e.g. surface brightness profile 
recevied from Hubble Space Telescope, then we can find the more precisely 
dependence of anisotropy parameter on radial distance in the concrete observed 
globular cluster. Further we plan to study nostationary effects. 

References 

Agekian, T. A., k Petrovskaya, I. V. 1962, Uchen.zap.LGU, 307, 187 
Bettwieser, E., & Spurzem, R. 1986, A&A, 161, 102 
Harris, W. E., & van den Berg, S. 1984, AJ, 89, 1816 

King, I. R., Hedemann, E. Jr., Hodge, S.M., & White, R. E. 1968, AJ, 73, 456 

Louis, P. D. 1990, MNRAS, 244, 478 

Meylan, G., & Heggie, D. C. 1997, A&ARew, 8, 31 

Peykov, Z. I., & Roussev, R. M. 1986, AZh, 63, 483 

Peykov, Z. I., & Roussev, R. M. 1988, AZh, 65, 317 

Solanes, J. M., & Salvador-Sole, E. 1990, A&A, 234, 93 

Spurzem, R. 1996, in IAU Symp. 174, Dynamical Evolution of Star Clusters - 
Confrontation of Theory and Observations, eds. P.Hut k, J.Makino (Dordrecht: 
KAP), 111 

Takahashi, K. 1996, in IAU Symp. 174, Dynamical Evolution of Star Clusters - 
Confrontation of Theory and Observations, eds. P.Hut Sz J.Makino (Dordrecht: 
KAP), 91 



